Astron. Nachr. / AN VVV, No. ii, 000-|6](yyyy) / DOI This.is/not.aDOI 



Current methods for analyzing light curves of solar-like stars 

J. Ballot* 

Laboratoire d'Astrophysique de Toulouse-Tarbes, Universite de Toulouse, CNRS, 14 av. E. Belin, 31400 Toulouse, France 

Received 1 May 2010, accepted 30 June 2010 
Published online later 



o 

(N 
< 

m 

6 

c3 



in 

On 

q 

o 
o 



X 



Key words stars: oscillations - stars: rotation - stars: activity - methods: data analysis 

CoRoT has allowed a quantitative leap for the solar-like-star seismology thanks to 5-month-long uninterrupted timeseries 
of high-precision photometric data. Kepler is also starting to deliver similar data. Now, several F and G main-sequence 
stars have been analyzed. The techniques developed to interpret light curves directly inherit from the experience got on 
the Sun with helioseismology. I describe in this review the methods currently used to analyze these light curves. First, 
these data provide an accurate determination of the stellar rotation rate. This is possible thanks to the magnetic activity of 
stars. The power spectra of light curves put also constraints on the stellar granulation, which can be directly compared to 
3-D stellar atmosphere models; this shows still unexplained discrepancies. I then detailed a standard method for extracting 
p-mode characteristics (frequency, amplitude and lifetime). CoRoT has revealed unexpected short life times for F stars. 
Last, I also discuss errors and biases of mode frequencies, especially the ones due to the simplified description of the 
rotation generally used. 
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1 Introduction ties, space-based missions perform high-precision photom- 

etry measurements. Acoustic modes generate small fluctua- 
In solar-like stars, acoustic (p) modes can be excited whereas| tions in the light curve which are analyzed. 

Thanks to their long duration, their high duty cycle, and 
their good signal-to-noise ratio, data sets obtained by space 
missions can be analyzed with techniques inherited from 
helioseismology. In this paper, we focus on the classical 
way to analyze such long photometry data sets for solar-like 
stars. We will illustrate these standard analysis techniques 
with recent results obtained with CoRoT. 



they are normally stable in regards to classical processes 
such as the well-known k mechanism. Damped modes are 
stochastically excited by the turbulent motions in their con- 
vective envelop. The resulting amplitudes are low. Exclud- 
ing the Sun, the first claim of a detection of solar-like p- 
mode oscillations ha s concerned Procyon . A power excess 
ha s been report e d by | Brown et al.1 d 199 lb . lately confirmed 



bv lMartic et al.l i 19991) . Nevertheless, the asteroseismology 
of solar-like st ars has really started after the o bservations 
of a Cen A by iBouchv & Carrier! d200ll l2002h which has 
revealed a very clear comb structure in the oscillation spec- 
trum. About 20 low-degree p modes in the spectrum have 
been identified. After the Sun, it was the first opportunity to 
do seismological stu dies and to test models for this kind of 
stars (see for instance 



2 Non-seismic information in seismic data 



2005; Thevenin et al 



High-precision photometry timeseries reveal more than acous-| 
tic modes in solar-type stars. Indeed, the luminosity fluctu- 
ates not only because of intrinsic oscillations, but also be- 
Eggenberger et alJ2004lMiglio & Montj 9Mnl of magnetic activity and surface motions. We discuss 
2002' Thoul et al [2003 Yildizl2007l) in this section what is the information we can obtained about 

the rotation of stars and their surface convection. 



Today, several F-, G-, K-type dwarf and sub-giant stars 
have revealed p-mode oscillations thanks to accurate radial- 
velocity measurements performed with high-resolution spec-| 2.1 Rotation 
trograph. To avoid the day-night alternation which gener- 

ates gaps in data and then windowing effects, one can run The first visible expression of the magnetic activity of the 
multi-site campaigns, such as those ded icated to Procyon Sun 18 the Presence of spots visible on the photosphere. 
dArentoft et al.l2008HBedding et alfeOld) . or use space-based| S P ots are 8 ood tackers " and historically the first ones - 
instruments which can follow continuously targets during of the solar rotation. If we use this solar description as a 
several weeks, several months, or even more. P modes in Paradigm, such spots should be present on the surface of 
solar-like star s have then been ac curately measured by CoRoT| solar - like stars with convective envelop which develop a 
dMicheletalibOOi and Kepler dChaplin et alfcoich . Con- magnetism similar to the Sun. The presence of spots, but 

also of bright areas called plages, creates modulations in 
the light curve with a period equal to the rotation period, 
e-mail: jballot@ast.obs-mip.fr Indeed, a magnetic structure creates such a modulation if 



trary to ground-based observations measuring radial veloci- 
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Fig. 1 Power spectrum in a low frequency range of Fig. 2 Power spectrum of a CoRoT target in gray. In 



the light curves of tw o CoRoT targets HD 181420 (af- 
ter iBarban et al.ll2009h and HD 181906 (after IGarcfa et al 



120091) showing the rotation frequency. In gray, the spectrum 
is oversampled with a factor of 4. 



black, the same plot smoothed over 100 points. Red dot- 
ted lines show the three components of the background (see 
text) fitted on the data, the red dashed line is the sum. The 
solid red line includes the p-mode contribution modeled 
with a Gaussian function. 



it is alternately visible or hidden on the far side. Due to 
the inclination i of the stellar rotation axis, structures too 
close to the poles are always visible or hidden and do not 
produce a modulation. Moreover, the lifetime of structures 
must be long enough, compared to the rotation period, to 
create repetitive patterns. Such modulations can be found in 
the power spectrum of stellar light curves. Contrary to rota- 
tion velocities measured in spectroscopy with line broaden- 
ing, these direct measurements of rotation periods is decou- 
pled from sinz. 

Figure[T]shows the power spectrum in the low frequency 
range for two CoRoT targets. Clear peaks appear as the sig- 
nature of the rotation, bu t patterns are diffe rent for the two 
stars. For HD 181420 dBarban et al.ll2009l) . the patt ern is 
broad with numerous peaks, whereas for HD 181 906 (IGarcfa 
2009), the power is concentrated in two peaks. We inter- 
pret the broadening as an effect of differential rotation in 
latitude. A spot tracks the rotation at the latitude it has ap- 
peared. When there are spots at various latitudes and a no- 
ticeable differential rotation, we observe a broad pattern in 
the power spectrum. On the other hand, a narrow peak can 
be due either to a weak differential rotation rate or to a small 
number of spots. For deeper analyses, it is also possible to 
directly fit the light curve with spot mo dels (see for instance 
Lanza et al.ll2006llMosser et alj|2009al) . 



It is worth noting that these observations give an in- 
formation on the presence of magnetic structures, even on 
small scales. It is then possible to have an indication of the 
stellar magnetic activity, according to the coverage of the 
surface with spots or plages. This information is comple- 
mentary to observations of large-scale m agnetic fields with 
spectropolarimetry (e.g. lPetit et alj|2008 ). 



2.2 Power spectrum of convection and background 

When we analyze the power spectrum of solar-like stars 
at higher frequency, we identify several components illus- 
trated in Fig. [2] The p-mode contribution is visible in this 
example as a bump around 2000 /iHz. The rest forms the 
background we decompose classically in three components 
as following: 

1 . a white noise, mainly due to photon noise, dominating 
the spectrum at the highest frequencies; 

2. a low-frequency component, which can be modelled with| 
a power law, resulting to activity effects and slowly vary-| 
ing trends; 

3. the contribution of the surface con vection, i.e. the gran- 
ulation, described with a model of lHarveyl (119851) : 

B c {u)= Jf l , (1) 

where v is the frequency, and a c , r c and a c are 3 free 
parameters. r c represents the characteristic time of the 
granulation. 

some authors (see iKarof j|201ob add also another contribu- 
tion due to faculae. Faculae are brights points visible at the 
solar surface that should also exist on other solar-type stars. 

For more than 10 years, the granulation patterns ob- 
served at the surface of the Sun are correctly reproduced by 
3-D models of the solar a tmosphere performed w ith radia- 
tive hydrodynamics code (IStein & Nordlundll998l) . Similar 
models have been p erformed for other stars. For instance, 



Ludwig et alJ (|2009) have performed models for the CoRoT 



target HD 49933 and they have compared the power spec- 
trum deduced from their simulations to the CoRoT obser- 
vations. When the granulation timescale is correctly repro- 
duced, the amplitude of the convective background is sig- 
nificantly different within a factor of 2-3 for this star. This 
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discrepancy has still to be explained. All this new bench 
of data offers an interesting opportunity to constrain this 
kind of 3-D models, which have already proved their values 
with their capability to consistently reprodu ce solar absorp- 
tion line profiles (e.g. lAsplund et al.ll2000l and following 
works). 

3 Mean seismic information 

For faint stars, due to the low signal-to-noise ratio, p modes 
are generally seen as a small power excess over the back- 
ground, but not as identifiable peaks. Nevertheless, it is still 
possible to extract seismic information. The most obvious 
information are v mstX , A max , and Av, which denote respec- 
tively the frequencies and amplitudes of the most powerful 
modes, located at the center of the power excess bump, and 
the large frequency separation (see fo r instance the CoRoT 
target HP 175726 a nd HD 170987 in IMathur et ai]l2010at 
Mosser et al.l2009 b). However, the estimates of A max are to 
be cautiously considered. Due t o a possible overestim ation 
of the background contribution. IMathur et a l. (2010b) have 
shown that, using their standard method, A max is systemat- 
ically underestimated. Detecting a power excess is not suf- 
ficient: it has to be shown that it is produced by the p modes 
and does not originate from another phenomena such as fac- 
ulae, as ment i oned in Sect. 12. 21 This problem is discussed in 



the frequency resolution sufficiently fine, the autocorrela 
tion function recovers the oscillatio n in Avjv) produced by 



the se cond helium ionisation (e.g. iMonteiro & Thompson 



2005). Such an information provides a valuable quantity to 



Mathur et al. ( 2010al) . To have a solid clue of its seismolog- 



ical origin, we must detect a periodic spacing in this power 
excess. Such a period spacing proves the presence of a comb 
structure hidden in the noise. By performing the autocorre- 
lation -or the spectrum- of the power spectrum restricted 
in a range around v max , we then determine the mean large 
separation {Av). 

In the past few years, several pipelines have been devel- 
oped to ex tract these global seismic parameters in an auto- 



constrain models more precisely. 

4 Individual p-mode determination 
4.1 Modeling the spectrum 

In the best cases, when the signal-to-noise ratio is suffi- 
ciently high, p modes are individually visible in the spec- 
trum. Acoustic modes in solar-like stars are considered as 
damped modes continuously and stochastically re-ex cited 
by the near-surface convection (e.g. Ku mar et al . 1988). The 
power spectrum of such a mode is a Lorentzian profile mul- 
tiplied by a stochastic noise following a exponential distri- 
bution. This approach is valid for the Sun. Although asym- 
metric Lor entzian profiles are gene rally considered in the 
solar case (lAbrams & Kumarlll996l) . nowadays, we totally 
neglect profile asymmetries for other stars due to the coarser 
frequency resolution and the lower signal-to-noise ratio. How-| 
ever, with longer observations, such as the Kepler long runs, 
the asymmetry will become a relevant parameter to fit, and 
will be needed to avoid biases in frequency determination. 
We also assume a low rotation rate il, without noticeable 
differential rotation. The rotation is treated as a first-order 
perturbation, and the Coriolis force effect is neglected since 
the mode frequencies ar e large compared to the rotation 
frequency (lLedouxlll95ll) . This assumption is discussed in 
Sect. [5] Thus, the frequency of a mode with radial order n, 
degree I and azimuthal order m is v n ,i,m = Vn,l — Jny 5 , 
where v n j is the frequency in absence of rotation and v s m 
v T0 t = S7/27T is the rotational splitting. 

Thus, we model the power spectrum S(v) of a light 



*" >■"- "«"-- h™"~+ r^^jw-w*^, y ^ curve as following 

matic way ( Hekker et al. 2010: IHuber et al.l2009tlKallinger et al. ] „, , _ > p , v» .,, » 
2010l: lMathur et alfeOlObh ."Automatic procedures are needed] ^ _ [ [y) ^ V >> X [y) 



today to efficiently exploit the observations of big sets of 
stars as those produced by Kepler. 

The determination of the mean large separation (Av), 
coupled with effective temperature measurements, provides 
an estimate of the stellar radius within an ac curacy of a few 
percents as discussed in lStello et al.l (12009). Of course, we 
then assume that stellar models are perfectly known. Such 
constraints do not improve our view of stellar physics nor 
refine our stellar evolution models, but are really useful for 
instance to characterize planet transits. 

To go beyond the mean large spacing {Av), the auto- 
correlation function (ACF) with a narrow filter is a valuable 
tool. The ACF can be computed as the Fourier spectrum 
of the filtered oscillation spectrum (iRoxburgh & Vorontsov 



(2) 

where B is the background described in Sect 12. 21 J\f a ran- 
dom function following a exponential law, and 



1 



P(v) 



H n ,iai, m {i) 



y y 

^ 1 + 2(i/ - v n .i + mv s ) Y n> {\ 

n.l m— — I 



(3) 



H n ,i are the modes heights, r n ; their widths - linked to the 
lifetimes t„.z = (7rr n! i) _1 - and ai, m (i) is the amplitude 
ratio, expressed as a function of the inclination angle: 



ai, m (i) 



\l 



-CFT(cosi)) 5 



(4) 



2006). If the bandpass filter applied to the oscillation spec- 



where P™ are the associated Legendre functions (e.g. Gizon & Sola nki] 
120031) . Moreover, we generally assume that the width is a 
slow function of the only frequency T n> i = T(v n .i). Simi- 
larly, we assume for the height, 

H n ,i = {Vi/V ) 2 H (v n .j), (5) 



wh ere Vj is the mode visibility (f or a generally discussion, 
see iToutain & Gouttebroze 1993). In other words, the en- 



trum is narrow enough, by scanning the p-mode frequency 
range, it is possible to recover variations of Av with fre- 
quency (for more detai ls on this technique, see lMosser & Appom nghaaori|ipartition between modes of close frequency is as- 
2009t [Roxburgh ^009j). If the noise is sufficiently low and sumed. Vi is rapidly decreasing with I and becomes very 
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small for I > 2. Thus, in practice, only modes I = 0, 1 
and 2 (and, sometimes, I = 3) are considered. Under very 
general assumption, it is worth noting we have the relation 



k 



V, 



l=2k 



k 



v, 



;=2fe+i 



(6) 



A demonstration is proposed in appendix. 

To recover the values of the free parameters, the model 
is fitted to the data using a maximum likelihood estima- 
tor (MUE} 2 _as_J^Jias been done for the Sun (see method 



Appourchaux et alj 1 1998). Nevertheless, it is now very 



common to consider Bayesian priors on the parameters and 
then to use maxim um a posteriori estimators (MAP, e.g. 
Gaulm e et al. I l2009h . In Bayesian approaches, Monte Carlo 
Markov chain s (MCMC) have sh own their usefulness on 



o/=i 
ni=2 



from PDF 
Kramer-Rao 
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CoRoT data dBenomar et al. I l2009l) : when classical estima- Fig. 3 Errors of mode frequencies v n j for HD 49933 after 



to the real errors determined on the PDF by Ben omar et al. 



tors only return the value of a maximum of probability (which| 180 days of observations. Black lines/symb ols correspond 
can be only local in the large parameter space to explore), 
MCMC can return the marginal probability density func- 
tion (PDF) for each parameter. PDFs sometimes appear to 
be non-normal, or even multimodal. 



4.2 Errors on parameter estimation 

The best way to determine the errors of parameters is to 
derive them from the PDF computed with MCMC. When 
one use s a MLE, the errors are es timated from the Hessian 
matrix dAppourchaux et al.ll 1998b . These estimations of er- 
rors are based on the Kramer-Rao theorem. Thus we have 
to keep in mind the following important points: 

1 . it is only a lower limit of the error; 

2. it is only asymptotically valid, so the PDF is supposed 
to be normal (that is not always the case with several 
mont hs of data, especially for i and v s , 



Ballot et al 



2008); 



3. the Hessian matrix should be computed for the real val- 
ues of the parameters, and not for the estimated values 
as done during a real fit. 

Some general properties of errors have been st udied last 
years. Concerning the inclination and the splitting, Ball ot et al 



(2006) have shown the strong correlation between the de- 
termination of i and v s , when v s < T (this is generally the 
case). Thus the errors of i and v s become huge, but pro- 
jected splittings v* = v s sini are better determined. 

The Kramer-Rao limit for v n x has be determin ed first 
for a single Lorentzian profile by iLibbrecht ( 1992b . It has 
been generaliz ed for multiplets in the solar configuration 
(i = 90°) by Toutain & Appourchaux! (1 1994b and for any 
inclination by lBallot et al.1 (120081) . The errors of v n j are: 



1 r r 



-fl(0n,hX s ,i), 



(7) 



4tt T 

with T the observation duration, = B(v n j)/ H n i the 
noise-to-signal ratio, and x s = 2v s /T n i the reduced split- 
ting. The functions /; are described in iBallot et alj (2008) 
and varies typically between 0.5 and 1.5 for C0R0T stan- 
dard targets. 



d2009h . Red lines/symbols show Kramer-Rao limits deter- 
mined with Eq. [7] 

4.3 The F-type star HD 49933 

The first solar-like star ob served by C0R0T was HD 49933 
dAppourchaux et al.l 120081) . It has revealed modes broader 
than expected. Mode widths are larger than the small sep- 
aration 0-2, thus / = 2 modes are blended with I = 
modes. Moreover, due to Eq. [6] the power of I = 1 modes 
is very close to the sum of the power of nearby I = 
and / = 2 modes. The modes are then not easily identi- 
fie d. With 60 days of data, th e mode identification favored 
by lAppourchaux et al.l ( 20081) has appeared to be wrong. It 
has been possible to undoubtedl y distinguish even and odd 
modes with longer observations ( Benomar et alj |2009). 

Such high mode widths has consequences on frequency 
accuracies, according to Eq.|7] Moreover, errors are increased^ 
by the blending of I = and 2 modes. Indeed, Eq.|7]is de- 
rived for isolated modes; when two modes are blended, it 
appears correlations in the determination of the two mode 
frequencies, that tends to increase the errors. This is illus- 
trated on Fig. [3] For HD 49933, the real errors of I = 1 mode 
frequencies, determined from the PDF, are almost equal to 
the Kramer-Rao limit given by Eq. |7](^ 0.5 ^Hz). How- 
ever, the errors for even modes, especially for I = 2 modes, 
are increased. Accurate determinations of I = 2 modes for 
this kind of F stars appear to be very difficult, and require 
year-long observations, as Kepler will do. 

4.4 Mode lifetimes 

Mode lifetimes appear to decrease when effective temper- 
atures T c ff increase. Using theoretical models for damping 
processes and some sei smic observat i ons b ased on radial 
velocity measurements, Chaplin et al. (2009) have found a 
strong depen dence r oc Tip- Rece ntly, by considering CoRoTj 
observations, Bau din et al. (2010) have even found a stronger] 
dependence V oc T^. This is different from what they have 
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Table 1 Perturbative coefficients representative of a 
HD 49933-like star used in this paper. C is a natural number, 
D 1 and D 2 are expressed in /iHz" 1 . These coefficients have 
been estimated by using results from Reese et al. ( 20061) and 
correspond to n ~ 10. 
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found for solar-like oscillations in red giants where T is 
weakly sensitive to T e s. 

Such a strong dependence on T c s has still to be ex- 
plained. Progresses should be done to understand the damp- 
ing processes. Nowadays, there is not cle ar consensus and 
differe nt mechanisms have been proposed. Goldreich & Ku marJ 
(119911) have shown that radiative losses and turbul ent vis- 
cosity pl ay a major role in mode damping, whereas Gough 




dl980b or lBalmforthl (1 1992b found that the damping is dom- 
inated by the modulation of turbulent pressure. More re- 
cently, iDupret et al. (2006) suggested the dominant contri- 
bution is the perturbation of the convective heat flux. 



5 Effects of rotation on the oscillations 

Perturbative treatments are valid only for low rotation rates. 
Is it the case in practice? As an illustration, we consider 
HD 49933. Its rotation frequency is v mt 3.4 /iHz, i.e. 
n w 0.05^ with n K = y/GM/R 3 the Keplerian fre- 
quency, derived fr om estimated mass M and radius R of the 
star. According to lReese et al. (2006), in this regime, pertur- 
bative treatments of the rotation are still valid, but second 
order terms must be considered. Frequencies are then ex- 
pressed 

v n ,i, m = v n j + {l-C n j)mv rot + {D] ll +m 2 D 2 nl )i> 2 ot {W^ 
according to ISaiol dl98lh . I have crudely estimated pertur- 
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Fig. 4 Synthetic spectra (black lines) showing three con- 
secutive modes I = 2, and 1, for HD 49933-like stars with 
different inclination angles i — 20, 60 & 90° (top to bottom 
panels). Mode frequencies are computed with Eq.|8] Green 
ticks indicate the frequencies without rotation, black ticks 
the frequencies with rotation (solid, dotted and dashed ticks 
correspond to \m\ = 0, 1 & 2 components). Red lines show 
the fit of Eq.[3]and red ticks indicate the fitted frequencies. 



bative coefficients C n ,i, -D 1 z and D 2 , from a poly tropic 



model by using results from iReese et al.l (2006) (see Ta- 
ble Q~|i. As expected, the Coriolis correction C n .i is very 
small and can be neglected. However the second-order terms| 
modify the frequencies by several /iHz. Multiple ts are shifted,| 
and they are no more symmetric relative to the m = 
component. Moreover, this estimation is optimistic, since 
the values used here correspond to low-order modes (n s» 
10): for higher orders, second-order terms should even be 
greater. 

Since this asymmetry is not taken into account in the 
fitted spectrum model (Eq. [3), the fitted parameters are bi- 
ased. To estimate these biases, I have computed synthetic 
spectra using the 2nd-order perturbative formulation (Eq.[S]l 
and performed classical fits. Results are shown in Fig. |4] 

Biases for frequencies v n j vary according to the incli- 
nation, typically from ~ 1.3 to 2.6 /iHz. For low angles 
(i < 30°), fits are dominated by the m = components, 



so the fitted frequencies are v, 



(fit) 



Dl 



n,i"rot 



. For 



high angles (i > 80°), spectra are dominated by the sec- 
toral modes \m\ = I, thus the fitted frequencies are 



v n ,i + (D^ i + l 2 D 2 n ] )v 2 ot . Fitted rotational splittings v s 
deviate also from the value of v mt , especially at low an- 
gles: the bias is around 0.3 /iHz for our case with i = 20°. 
However, for i > 80°, the bias vanishes: the sectoral modes 
dominate and they are still separated by 2lv mt . Concerning 
the other parameters (i, H and T), the biases are small and 
negligible relatively to the errors. 

Biases are significant for frequencies and splittings in 
stars like HD 49933. Second-order terms must then be con- 
sidered. Ideally, it should be included in the fitted mod- 
els, but it increases the number of free parameters and risk 
destabilising the fit. Nevertheless it must be taken into ac- 
count for the interpretation of frequencies. 
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6 Conclusion 

Last years, seismic data of unparalleled quality have been 
produced for solar-like stars, thanks to CoRoT, and now, 
thanks to Kepler. The delivered high-precision photomet- 
ric timeseries put a lot of news constraints on stars. First, 
we get accurate measurements of rotation rates and good 
indications for magnetic activities. Next, the background of 
the spectrum as well as mode amplitudes and lifetimes con- 
strain the turbulent dynamics of the near-surface layers, and 
it offers new capabilities to test 3-D stellar atmosphere mod- 
els. The extracted frequencies must now be interpreted in 
terms of constraints on the stellar structure, especially to 
improve the treatment of mixing processes. Nevertheless, 
the rotation has to be carefully considered, not only because 
of the effects induced on the structure and the evolution of 
stars, but also because frequencies are significantly affected 
even for typical F main-sequence stars. 

Acknowledgements. The author acknowledges the support of the 
Agence National de la Recherche through the SIROCO project. 
This work was supported by the European Helio- and Asteroseis- 
mology Network (HELAS), a major international collaboration 
funded by the European Commission's Sixth Framework Programn 

References 

Abrams, D., Kumar, P.: 1996, ApJ 472, 882 
Appourchaux, T, Gizon, L., Rabello-Soares, M.: 1998, A&A 132, 
107 

Appourchaux, T, Michel, E., Auvergne, M., et al.: 2008, 
A&A 488, 705 

Arentoft, T, Kjeldsen, H., Bedding, T. R., et al.: 2008, ApJ 687, 
1180 

Asplund, M, Nordlund, A., Trampedach, R., Allende Prieto, C, 

Stein, R. R: 2000, A&A 359, 729 
Ballot, J., Appourchaux, T, Toutain, T, Guittet, M.: 2008, 

A&A 486, 867 

Ballot, J., Garcia, R. A., Lambert, P.: 2006, MNRAS 369, 1281 
Balmforth, N. J.: 1992, MNRAS 255, 639 
Barban, C, Deheuvels, S., Baudin, F, et al.: 2009, A&A 506, 51 
Baudin, F, Barban, C, Belkacem, K, et al.: 2010, A&A, submit- 
ted 

Bedding, T. R., Kjeldsen, H., Campante, T. L., et al.: 2010, 
ApJ 713, 935 

Benomar, O., Baudin, F, Campante, T. L., et al.: 2009, A&A 507, 
L13 

Bouchy, F, Carrier, F: 2001, A&A 374, L5 
Bouchy, F, Carrier, F: 2002, A&A 390, 205 
Brown, T. M., Gilliland, R. L., Noyes, R. W., Ramsey, L. W.: 1991, 
ApJ 368, 599 

Chaplin, W. J., Appourchaux, T, Elsworth, Y., et al.: 2010, 
ApJ 713, L169 

Chaplin, W. J., Houdek, G, Karoff, C, Elsworth, Y., New, R.: 

2009, A&A 500, L21 
Dupret, M. A., Barban, C, Goupil, M., et al.: 2006, ESA-SP 624, 

97 

Eggenberger, P., Charbonnel, C, Talon, S., et al.: 2004, A&A 417, 
235 

Garcia, R. A., Regulo, C, Samadi, R., et al.: 2009, A&A 506, 41 



www.an-journal.org 



Gaulme, P., Appourchaux, T, Boumier, P.: 2009, A&A 506, 7 
Gizon, L., Solanki, S. K.: 2003, ApJ 589, 1009 
Goldreich, P., Kumar, P.: 1991, ApJ 374, 366 
Gough, D.: 1980, LNP 125, 273 
Harvey, J.: 1985, ESA-SP 235, 199 

Hekker, S„ Broomhall, A.-M., Chaplin, W. J., et al.: 2010, MN- 
RAS 402,2049 

Huber, D., Stello, D., Bedding, T. R., Chaplin, W. J., Arentoft, T, 
Quirion, P.-O., & Kjeldsen, H.: 2009, CoAst 160, 74 

Kallinger, T, Weiss, W. W., Barban, C, et al.: 2010, A&A 509, 
A77 

Karoff, C: 2010, AN, this volume 

Kumar, P., Franklin, J„ Goldreich, P.: 1988, ApJ 328, 879 
Lanza, A. F, Messina, S., Pagano, I., Rodono, M.: 2006, AN 327, 
21 

Ledoux, P.: 1951, ApJ 114, 373 
Libbrecht, K. G.: 1992, ApJ 387, 712 

Ludwig, H., Samadi, R., Steffen, M., et al.: 2009, A&A 506, 167 
Martic, M., Schmitt, J., Lebrun, J., et al.: 1999, A&A 351, 993 
Mathur, S., Garcia, R. A., Catala, C, et al.: 2010a, A&A in press 
Mathur, S., Garcia, R. A., Regulo, C, et al.: 2010b, A&A 511, 
A46 

Michel, E., Baglin, A., Auvergne, M., et al.: 2008, Science 322, 
558 

Miglio, A., Montalban, J.: 2005, A&A 441, 615 
Monteiro, M. J. P. F. G., Thompson, M. J.: 2005, MNRAS 361, 
1187 

Mosser, B., Appourchaux, T: 2009, A&A 508, 877 
Mosser, B., Baudin, F, Lanza, A. F, et al.: 2009a, A&A 506, 245 
Mosser, B., Michel, E., Appourchaux, T, et al.: 2009b, A&A 506, 
33 

Petit, P., Dintrans, B„ Solanki, S. K, et al.: 2008, MNRAS 388, 
80 

Reese, D., Lignieres, F, Rieutord, M.: 2006, A&A 455, 621 

Roxburgh, I. W.: 2009, A&A 506, 435 

Roxburgh, I. W., Vorontsov, S. V.: 2006, MNRAS 369, 1491 

Saio, H.: 1981, ApJ 244,299 

Stein, R. F, Nordlund, A.: 1998, ApJ 499, 914 

Stello, D., Chaplin, W. J., Bruntt, H., et al.: 2009, ApJ 700, 1589 

Thevenin, F, Provost, J., Morel, P., et al.: 2002, A&A 392, L9 

Thoul, A., Scuflaire, R., Noels, A., et al.: 2003, A&A 402, 293 

Toutain, T. & Appourchaux, T. 1994, A&A, 289, 649 

Toutain, T, Gouttebroze, P.: 1993, A&A 268, 309 

Yildiz, M.: 2007, MNRAS 374, 1264 



A Odd and even mode visibilities 

The visibility Vj is given by the relation dGizon & Solankill2003l) : 

Vi = 2n [ Yfi^Wirindfi (Al) 
Jo 

where F ; m denotes the spherical harmonics and W is a function 
indicating the contribution to the total light flux of an element of 
the stellar disk; we have assumed that W depends only on the dis- 
tance to the limb fi. We rewrite the equation 

Vi = v/7r(2i + 1) / Pi(n)f(vW (A2) 
Jo 

with f(fi) = W(fi)fi and Pi the Legendre polynomials. We define 
/ e (fi) the even function over the interval [—1,1] such that /„ (/j) = 
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/(/i) Wfi G [0, 1] and f (fj) the odd function over the interval 
[-1, 1] such that f a (fi) = f(p) V M G [0, 1]. 

Legendre polynomials form a orthogonal basis of the Hilbert 
space Ij2([— 1, 1]) where is defined a dot product 

(f,g) = f f(x)g(x)dx. (A3) 

Legendre polynomials verify the relation 

As f e and f are obviously elements of L2([— 1, 1]), we decom- 
pose them on the Legendre basis and take into account their parity: 

oo 

f e (x) = ^a 2k P 2 k{x) (A5) 

fc=0 
oo 

fo(x) = Y.^+iP2 k +i{x) (A6) 

fc=0 

Thus for I = 2k and I = 2k + 1, we obtain respectively 



V 2k = ^(4* + 1) i P 2k {ti)f e {ti)dti 



(AT) 



" "2fc (A8) 



4fe + 1 



1 

VWi = V7r(4fe + 3)- y P 2 k+i(n)fe(n)dn (A9) 
T /Wi (A10) 



4fc + 3 

We consider now (/ e , / e ) the square of the norm of / e . On one 
hand, we get by using its decomposition on the Legendre basis 

(/.»/•> =E^*7TT 7- (A11) 



'4fc + 1 

fc=0 

On the other hand, 

\2 



•/ 



(/«,/«)=/ U(xydx = 2 f{xYdx (A12) 



then 



o° 2 /-l 

fc=0 • / ° 



By considering (/„, /„), we show similarly that 

4fc + 3 



oo o2 /-l 

E^7=/ /(^- (A14) 



fc=0 

Hence we have 

,2 °° fl2 



^4fc + l &4k + 3 J 
Finally, by injecting Eqs |A8l & |A10| in Eq. |A15l we get 

oo oo „ ^ 

E V ™ = E V »+i =lr WOO W (A16) 

I — n n ■'0 
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